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Titanate compounds have been recognized as key materials for understanding the coupling of 
magnetism and orbitals in strongly correlated electron systems. In the perovskite Ti oxide _RTiC>3 
(R =trivalent rare-earth ions), which is a typical Mott-Hubbard insulator, the Ti t^g orbitals and 
spins in the 3d 1 state couple each other through the strong electron correlations, resulting in a rich 
variety of orbital-spin phases. One way of controlling the coupling is to change the tiltings of the 
TiC>6 octahedra (namely the GdFe03-type distortion) by varying the R ions, through which the 
relative ratio of the electron bandwidth to the Coulomb interaction is controlled. With this control, 
these Mott insulators exhibit an antiferromagnetic-to-ferromagnetic (AFM-FM) phase transition, 
which has turned out to be a consequence of rich orbital physics in these materials. The origin 
and nature of orbital-spin structures of these Mott insulators have been intensively studied both 
experimentally and theoretically. When the Mott insulators are doped with carriers, the titanates 
show touchstone properties of the filling controlled Mott transition. In this article, we first review the 
state of the art on the studies for understanding physics contained in the properties of the perovskite 
titanates. On the properties of the insulators, we focus on the following three topics: (1) the origin 
and nature of the ferromagnetism as well as the orbital ordering in the compounds with relatively 
small R ions such as GdTiOa and YTiOs, (2) the origin of the G-type antiferromagnetism and 
the orbital state in LaTiOs, and (3) the orbital-spin structures in other AFM(G) compounds with 
relatively large R ions (R =Ce, Pr, Nd and Sm). On the basis of these discussions, we discuss the 
whole phase diagram together with mechanisms of the magnetic phase transition. On the basis of the 
microscopic understanding of the orbital-spin states, we show that the Ti t2 g degeneracy is inherently 
lifted in the titanates, which allows the single-band descriptions of the ground-state and the low- 
energy electronic structures as a good starting point. Our analyses indicate that these compounds 
offer good touch-stone materials described by the single-band Hubbard model on the cubic lattice. 
From this insight, we also reanalyze the hole-doped titanates R\- x A x TiOz (A ^divalent alkaline- 
earth ions). Experimentally revealed filling-dependent and bandwidth-dependent properties and the 
critical behavior of the metal-insulator transitions are discussed in the light of theories based on the 
single-band Hubbard models. 
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I. INTRODUCTION 



After the discovery of high-T c superconductivity in the 
layered cuprate perovskites, strongly correlated electrons 
in the perovskite-type transition-metal oxides have at- 
tracted great interest 1] . Extensive studies have revealed 
astonishing transport phenomena upon carrier doping in 
these compounds. 

One of the important achievements of these studies 
is a discovery of the colossal magnetoresistance in the 
Mn perovskites. The manganites have the e g bands at 
the Fermi level similarly to the cuprate superconductors. 
However, it has turned out that the Mn 3c? orbitals play 
crucial roles on the electronic properties of the mangan- 
ites, which is in sharp contrast with the cuprate super- 
conductors in which the orbital degree of freedom is fully 
quenched because of both 3d 9 configuration of the mother 
compounds and large crystal-field splitting. While the 
strong electron correlations in the high-T c superconduc- 
tors have been studied based on several single-band mod- 
els including only the charge and spin degrees of freedom, 
those in the manganites have been studied by consider- 
ing the coupling of not only charges and spins but also 
orbitals, which revived a field of the condensed-matter 
physics, so called "orbital physics" 0. 



Perovskite Ti oxide is also one of the important ex- 
amples, and has attracted great interest. The titanates 
have the t^g bands at the Fermi level in contrast with the 
cuprates and manganites, and show very different prop- 
erties from these two compounds in spite of the struc- 
tural similarity. The perovskite titanate R\- x A x 'Y\0?, 
(i?=trivalent rare-earth ions, A=divalent alkaline-earth 
ions) shows a transition from insulator to metal upon 
carrier doping achieved by increasing x. Thus, this sys- 
tem has been focused on as a typical example of the 
filling-control metal-insulator transition system 
The metallic phase proved to be well described by the 
Fermi liquid picture in marked contrast with the super- 
conducting cuprates and the magnetoresistive mangan- 
ites. As the system approaches the metal-insulator tran- 
sition point from the metallic side, strong mass enhance- 
ment was experimentally found in the specific heat and 
the magnetic susceptibility. Anomalous filling-dependent 
behaviors near the transition point have been observed 
in various experiments, and have been discussed in con- 
nection to the nature of the Mott transition 0,Si- The 
origin of the very different behavior of the titanates from 
the cuprates and manganites has been a subject of very 
intensive studies, and the titanates have been recognized 
as an important class of materials in clarifying the strong 
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electron correlation effects. 

In order to clarify the nature of the metallic phase and 
the transition, it is also important to clarify the electronic 
structures of the end Mott-Hubbard insulators i?TiC>3. 
Similarly to the manganite case, the Ti 3c? orbitals have 
been recognized to be important for the electronic struc- 
tures in i?TiC>3 as well. In the Mott insulator, while 
the low-energy excitations of the charge degree of free- 
dom are quenched by the insulating gap, the orbital and 
spin degrees of freedom survive, and the couplings and 
interplays of these two due to the strong electron corre- 
lations give rise to various electronic structures and puz- 
zling phenomena. Coupling to lattice further enriches 
their interplays. 

Because of a variety of orbital-spin phases and their 
phase transitions due to these couplings and interplays, 
the perovskite titanate i?Ti03 has been considered as a 
key class of materials for the understanding of the or- 
bital physics. Because of the nominal three-fold degen- 
eracy of d 1 state under the cubic crystal field on Ti ions 
in the insulating titanate, it is expected that the com- 
plete understanding of the Ti perovskites contributes in 
clarifying roles of the orbital degree of freedom, and cou- 
plings between magnetism and orbitals in the strongly 
correlated electron systems. In addition, the knowledge 
about the couplings is expected to open a possibility of 
controlling transport, optical and magnetic properties by 
various external control parameters like magnetic fields, 
electric fields, light, pressure etc. 

In this article, we first discuss recent progress of stud- 
ies on the orbital-spin structures in the insulating ti- 
tanates i?TiO-3 and Lai-aYxTiOg. We point out most 
of the available experimental results are accounted for 
by careful analyses of the orbital-spin structures and ex- 
citations. In particular, the realistic crystal field lifts 
the nominal ti g orbital degeneracy, providing a firm ba- 
sis and justification of the description by a single-band 
Hubbard model. This indicates that the perovskite ti- 
tanates indeed offer touch-stone materials for the single- 
band Hubbard model. From this justification, we re- 
visit the bandwidth-control experiments for these com- 
pounds together with the theoretical analyses in detail 
and concludes that the single-band Mott transition is 
indeed realized in the titanates. In addition, we also 
reanalyze the experimental studies on the filling-control 
metal-insulator transition in the hole-doped systems by 
favorably comparing with theoretical results of single- 
band models. The limitation of the single-band descrip- 
tion in real titanates is also discussed. The rest of this 
article is organized as follows. We first introduce the ex- 
perimentally obtained magnetic phase diagram for the 
insulating titanates, and mention some puzzles existing 
in this phase diagram. In Sees 3-6, by reviewing a num- 
ber of attempts to solve these puzzles and the solutions 
provided by recent studies, we discuss the orbital-spin 
structures in the insulating titanates as well as the nature 
of their phase transition from a unified viewpoint. The 
bandwidth-control experiments and corresponding theo- 
retical studies are also discussed in Sec. 6. In Sec. 7, we 



FIG. 1: Crystal structure of i?Ti0 3 with the GdFe0 3 -type 
distortion. 

discuss the nature of the filling-control metal-insulator 
transition in the hole-doped systems based on both ex- 
perimental and theoretical results. Section 8 is devoted 
to the concluding remarks. 



II. MAGNETIC PHASE DIAGRAM FOR THE 
INSULATING TITANATES OTIOs 

First, we briefly review some fundamental properties of 
the perovskite Ti oxides. In the end compound i?Ti03, 
Ti 3+ has a t\ g configuration in which one of the threefold 
tig orbitals is occupied by an electron. 

The crystal structure of i?TiC"3 is a pseudocubic per- 
ovskite with an orthorhombic distortion (the GdFeC>3- 
type distortion) in which the TiC>6 octahedra forming 
the perovskite lattice tilt alternatingly. Note that with 
this distortion, the cubic TiOg octahedra themselves are 
hardly distorted in contrast with the Jahn- Teller distor- 
tion. In this structure, the unit cell contains four TiOg 
octahedra (sites 1-4) as shown in Fig. 2] The magni- 
tude of the distortion depends on the ionic radii of the R 
ions. With a small ionic radius of the R ion, the lattice 
structure is more distorted and the Ti-O-Ti bond angle is 
decreased more significantly from 180°. For example, in 
LaTiC>3, the bond angle is 157° (a&-plane) and 156° (c- 
axis), but 144° (afr-plane) and 140° (c-axis) in YTi0 3 [|. 
The ionic radii of La and Y ions are 1.17 A and 1.04 A, 
respectively. The magnitude of the distortion can also be 
controlled by using solid-solution systems. For example, 
by varying the Y concentration in Lai-^Y^TiOs, we can 
control the bond angle almost continuously from 157° 
(x=Q) to 140° (x=l) @]. 

It is widely recognized that the tilting of the TiO*6 oc- 
tahedra primarily controls the electron transfers between 
the neighboring Ti ti g orbitals mediated by the O 2p or- 
bitals, thereby reduces the ti g bandwidth with decreas- 
ing bond angle. Therefore, the GdFeO-3-type distortion 
has been recognized as a bandwidth-control mechanism, 
and the variations of the Mott-Hubbard gap and electron 
correlations in the insulating titanates have been exper- 
imentally studied by controlling the bandwidth through 
this distortion |8llall(]|. 

In addition to the bandwidth control, the tilting plays 
other important roles on the electronic structures. In 
fact, the GdFeO-3-type distortion controls the coupling 
and interplay of the orbital, spin and lattice degrees of 
freedom in the insulating titanates, and the system shows 
various orbital-spin phases and their phase transitions as 
we will discussed in the following. 

A magnetic phase diagram in the plane of temperature 
and the magnitude of the GdFeO-3-type distortion was 
experimentally obtained for i?Ti03, which exhibits an 
antifcrromagnetic (AFM)-to-ferromagnetic (FM) phase 
transition (see Fig. Supper panel) LaTiO-3 with 
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FIG. 2: Experimentally obtained magnetic phase diagram 
for i?Ti0 3 (upper panel, Ref. El IH), and that for 
Lai-zY^TiOa (lower panel, Ref. UU^) from Ref. In 
thes phase diagrams, Tn and Tc are plotted as functions of 
the ionic radius of the R ion and the unit cell volume, respec- 
tively. (The lines are guides for the eyes.) Note that the unit 
cell volume of Lai-^Y^TiOs is approximately proportional 
to the Y concentration (x), and well characterizes the mag- 
nitude of the GdFe03-type distortion |13|| . Here, the open 
symbols and the closed symbols indicate Tn and Tc, respec- 
tively. The values of the unit cell volume are deduced from the 
data obtained by the previous x-ray diffraction measurement 
(Ref. 0). 

the smallest distortion shows a G-type AFM [AFM(G)] 
ground state, in which spins are aligned antiferromagnet- 
ically in all x, y and z directions as shown in the inset of 
the upper panel of Fig. The magnetic moment is ^0.57 
Mb U3 ; which is reduced from the spin- 1/2 moment of 
1 fi B - The Neel temperature (T N ) is about 140-150 K. 
With increasing GdFe03-type distortion, Tn gradually 
decreases and is strongly depressed at SmTi03, subse- 
quently a FM ordering appears. In the significantly dis- 
torted compounds such as GdTi03 and YTi03, a FM 
ground state accompanied by a large Jahn- Teller distor- 
tion is realized. A similar phase diagram was also ob- 
tained for Lai^Y^TiOs (see Fig. flower panel) |slll3|. 

These phase diagrams provide us with several puzzling 
issues. First, strong suppressions of Tn and Tq around 
the AFM-FM phase transition point imply a continuous- 
type transition at T — 0. At first sight, this is a puzzling 
feature, because we expect the first-order transition be- 
tween completely different symmetry breakings at T = 0, 
and Tn and Tc should remain nonzero at the transition 
point. Clarification of the mechanism of this continuous- 
type transition has long been a issue of intensive studies. 

To clarify the nature and mechanism of this magnetic 
phase transition, it is essentially important to clarify the 
nature of the AFM(G) and FM phases. However, the ori- 
gin of the AFM(G) ordering in the compounds with weak 
GdFe03-type distortion has also been a puzzle. In the 
orbitally degenerate systems, there is a strong tendency 
towards a FM state accompanied by an antiferro-orbital 
ordering Q, while i?Ti0 3 with R =La, Ce, Pr, Nd, Sm 
show actually AFM(G) order at low temperatures. Here, 
"antiferro-orbital" order is defined as an order of orbital 
polarization with the staggered pattern in analogy with 
the staggered spin polarization in the AFM order. On 
the other hand, it turned out that the orbital state in 
the FM compounds with strong GdFe03-type distortion 
such as GdTiOs and YTi03 cannot be described by a 
simple antiferro-orbital ordering, which has also been a 
puzzle. 

As discussed above, there exist a number of puzzles 
in the experimentally obtained phase diagram for the in- 
sulating titanates. It has been widely recognized that 
to fully understand the electron correlation effects in 
the perovskite compounds, these puzzles have to be ex- 
plained by a convincing theory. In the following four sec- 



FIG. 3: Orbital ordering in YTi0 3 . (From Akimitsu et al., 
Ref. |H) 

tions (Sees. 3-6), we discuss the orbital-spin structures in 
the insulating titanates as well as the nature of the phase 
diagram. In Sec. 3, we discuss the orbital ordering and 
the origin of the ferromagnetism in the FM compounds 
such as GdTiOs and YTi03, which are located in the 
region of strong GdFe03-type distortion. In Sec. 4, the 
AFM(G) ordering and the orbital state in LaTiOs with 
the smallest GdFe03-type distortion are discussed. In 
Sec. 5, the orbital-spin structures in other AFM(G) com- 
pounds (R —Ce, Pr, Nd and Sm) in the weakly distorted 
region are discussed. Finally, we discuss the whole nature 
of the magnetic phase diagram and the phase transition 
in Sec. 6. 



III. ORBITAL ORDERING AND 
FERROMAGNETISM IN THE COMPOUNDS 
WITH STRONG GDFEO3-TYPE DISTORTION 

In this section, we discuss the orbital-spin structures 
in the FM compounds with relatively small R ions such 
as GdTi03 and YTi03. These compounds are located 
in the region of strong GdFe03-type distortion in the 
phase diagram of Fig. [21 Clarification of the origin and 
nature of the FM ordering is important for understanding 
the mechanism of puzzling second-order-like AFM-FM 
phase transition. Further, it is tightly connected with 
understanding the coupling between the magnetism and 
orbitals in the perovskite compounds. Hence, the issue 
has long attracted great interest, and a lot of studies 
have been done both theoretically and experimentally. 
We here overview the history of studies on this issue, 
and summarize recent progresses. 

The emergence of the FM ordering in the titanates 
cannot be explained by a single-band model, and it has 
turned out indispensable to take the 3d orbital degen- 
eracy into account to reproduce it. A model Hartree- 
Fock study based on the multiband d-p model succeeded 
in reproducing the FM ordering in YTi03, and showed 
that this FM state is accompanied by an orbital order- 
ing [fH ITU . The first-principles band calculations based 
on the generalized gradient approximation (GGA) and 
the local spin-density approximation (LSDA) also suc- 
ceeded in reproducing the FM state with the orbital or- 
dering in YTiOs [la . fl7| . In the predicted orbital order- 
ing, the wavefunctions of the occupied orbitals at sites 
1, 2, 3 and 4 are c\yz + c-ixy, c\zx + c%xy, c\yz — c^xy 
and c\zx — c-ixy (ci 2 + C2 2 = 1 and c\ ~ V0.5-V0.6, 
respectively (see also Fig. EJ - This orbital ordering was 
actually observed in YTiOs by several experiments such 
as NMR (Refs. Ijjjlll9l l2c | ) and polarized neutron scat- 
tering (Refs. l2ll l22l] ) . A resonant x-ray scattering exper- 
iment (Ref. 123) also succeeded in detecting this orbital 
ordering in YTi03 although we should mention that the 
mechanism of resonant x-ray scattering still has a con- 
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troversy El 111- 

In the orbitally degenerate systems, there is a strong 
tendency towards FM state with antiferro-orbital order- 
ing In the antiferro-orbital order- 
ing, the neighboring occupied orbitals are orthogonal to 
each other. This orbital-spin configuration is stabilized 
by the spin-orbital exchanges due to both Hund's-rule 
coupling and electron transfers. Indeed, the occupied or- 
bitals in YTiC>3 are nearly orthogonal in the afo-plane. 
However, the occupied orbitals along the c-axis are not 
necessarily orthogonal. In this sense, the emergence of 
the FM state in YTiC>3 cannot be understood straightfor- 
wardly as the ferromagnetism with antiferro-orbital or- 
dering. The origin of the FM ordering was not clarified 
in the above theoretical studies. Moreover, the puzzle of 
the AFM-FM phase transition was also not resolved in 
these studies. 

Recently, these issues were studied on the basis of an 
effective model in the insulating limit [33,|34(. This study 
showed that; 

• The observed orbital ordering is stabilized by the 
hybridizations between the neighboring e g and t 2g 
orbitals induced by the GdFe03-type distortion, 
which are prohibited without distortion. 

• The FM spin-exchange along the c-axis is realized 
by the super-exchange processes mediated by the e g 
orbitals due to the increased ti g -e g hybridizations 
in the region of strong GdFe03-type distortion. 

• In the FM phase, the FM exchange along the c-axis 
becomes weaker with decreasing GdFe03-type dis- 
tortion since the ti g -& g hybridization is decreased. 
Hence, a two-dimensional anisotropy appears near 
the AFM-FM phase boundary, which results in the 
decrease of Tq as the system approaches the AFM- 
FM phase boundary. 

Let us discuss some details of this study in the fol- 
lowing: First, we briefly refer to the model used. The 
starting point is the multiband d-p Hamiltonian which 
includes full degeneracies of the Ti 3d and O 2p orbitals. 
The Hamiltonian is given by 

H dp = Hdo + H p + H t d v + H tpp + ff on _ S ite, (1) 

The first and second terms denote the bare Ti 3d and 
the O 2p levels, respectively. The third and fourth 
terms represent the d-p and p-p transfers, respectively, 
which are determined by using the Slater-Koster pa- 
rameters [H |3(J. The values of the Slater-Koster pa- 
rameters were estimated in the cluster-model analysis of 
valence-band and transition-metal 2p core-level photoc- 
mission spectra [37t l38| and the first-principles band cal- 
culation [33. The effects of the lattice structure along 
with the lattice distortions are considered by modifying 
the transfer integrals. The fifth term expresses the on- 
site Coulomb interactions, which consists of the following 
four contributions; 

^on-site = H u + H u i + Hj + Hy . (2) 



FIG. 4: d-type Jahn- Teller distortion. 

FIG. 5: Level splitting of the 3d orbitals in the cubic crys- 
tal field and that in the d-type Jahn- Teller distortion. With 
the d-type Jahn- Teller distortion, the threefold degenerate t2 g 
level splits into twofold degenerate lower levels and a nonde- 
generate higher level at each Ti site. The twofold degenerate 
e g level also splits . The ways of the level splitting are different 
between sites 1, 3 and sites 2, 4. 

Here, H u and H u > denote the intra- and inter-orbital 
Coulomb interactions, respectively, and Hj and Hj> de- 
note the exchange interactions. The term Hj is the origin 
of the Hund's-rule coupling which favors the spin align- 
ment in the same direction on the same atoms. The term 
Hj' gives the fi-pair hopping between the 3d orbitals on 
the same Ti atom. These interactions are expressed by 
using Kanamori parameters, u, u', j and j', which satisfy 
the following relations; u = u' + 2j and j — j' [4(J, |4l| . 
Because of these relations, the Kanamori parametrization 
retains the rotational invariance. The values of these pa- 
rameters were also estimated by the cluster-model analy- 
sis of the photoemission spectra [37], ^ . By integrating 
out the O 2p orbitals in the path-integral formalism, we 
arrive at the multiband Hubbard model, which includes 
only Ti 3d orbitals and the transfer integrals between the 
Ti 3d orbitals mediated by the O 2p orbitals. By applying 
the second-order perturbational expansion with respect 
to the electron transfers between the Ti 3c? orbitals to 
thus obtained multiband Hubbard Hamiltonian, we de- 
rive an effective spin-pseudospin Hamiltonian H ST on the 
subspace of states only with singly-occupied t2 g orbitals 
at each Ti site in the limit of the strong Coulomb repul- 
sion. In this effective Hamiltonian, only spin and orbital 
degrees of freedom are incorporated, and the threefold t2 g 
orbitals are expressed by using the pseudospin-1 operator 
t. Further adding a crystal field term H cry , we obtain 
the following effective spin-pseudospin Hamiltonian; 

H e e. — H cry , + H ST . (3) 

In order to study the FM phase in the region of strong 
GdFe03-type distortion, we have to take into account the 
effect of d-type Jahn- Teller distortion, which was experi- 
mentally observed in the FM compounds such as GdTiOs 
and YTi03. With the d-type Jahn- Teller distortion, the 
TiOe octahedra are elongated along the y direction at 
sites 1 and 3, while at sites 2 and 4, they are elongated 
along the x direction (see Fig. . The Ti 3d orbitals di- 
rected along the elongated bonds are lowered in energy as 
compared with the orbitals along the shorter bonds be- 
cause of the reduction of the repulsive potential from the 
surrounding O ions. As a result, with the d-type Jahn- 
Teller distortion, the xy and yz orbitals are lowered at 
sites 1 and 3, and the xy and zx orbitals are lowered 
at sites 2 and 4 as shown in Fig. [S] The effect of the 
Jahn- Teller distortion can be treated in the multiband 
d-p Hamiltonian by tuning the Slater-Koster parameters 
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FIG. 6: Total energies of the AFM(A), FM and AFM(G) 
states near the AFM(A)-FM phase boundary are plotted as 
functions of the Ti-O-Ti bond angle, which are calculated by 
using the effective spin-pseudospin Hamiltonian in the insu- 
lating limit (see text). Inset: energy difference between the 
AFM(A) and FM solutions. (From Mochizuki and Iniada, 
Ref. H) 



FIG. 7: Second-order perturbational energy gains strongly 
depend on the spin configuration. Substantial transfer pro- 
cesses along the c-axis for the energy gains are illustrated for 
both parallel and antiparallel spin configurations. The angles 
#afm(A) and #fm2 are fixed at 45°. The cross symbol x rep- 
resents the forbidden transfer process. (From Mochizuki and 
Imada, Ref. Q) 



for the longer and shorter Ti-0 bonds because it is empir- 
ically known that the magnitude of the transfer integral is 
proportional to d~ 3 5 where d being the bond length [3(|. 

Total energies of various orbital-spin configurations 
were calculated by solving the self-consistent mean-field 
equations. The result showed that in the region of strong 
GdFe0 3 -type distortion, there exist AFM(A) and FM so- 
lutions accompanied by a certain type of orbital ordering 
as stable saddle-point solutions, and the AFM(A)-to-FM 
phase transition occurs with increasing GdFe03-type dis- 
tortion as shown in the inset of Fig. Note that the 
AFM(G) solution has much higher energies, and is un- 
stable relative to the AFM(A) and FM solutions in this 
region. 

As for the orbital state, the orbital ordering realized in 
the AFM(A) phase and that in the FM phase are similar 
to each other. We can specify the orbital states by using 
angles #afm(A) and 9fm as, 

site 1; cos8 x \xy) + sin 9 x \yz), 

site 2; cos8 x \xy) + sin9 x \zx) , 

site 3; — cos 8 X \ xy) + sin 9 x \yz), 

site 4; — cos 9 X \ xy) + sin 8 x \zx), 



(4) 



where x =AFM(A), FM. Both 6> A fm(A) and 9 FM take 
~ 45°, and the difference between them is very small. 
Namely, the orbital ordering in which -^(xy + yz), 

-j^{xy + zx), ^(-xy + yz) and -^(~xy + zx) are oc- 
cupied at sites 1, 2, 3, and 4, respectively, are realized in 
both AFM(A) and FM phases, and this orbital structure 
hardly changes through the AFM(A)-FM phase transi- 
tion. We refer to this orbital ordering as (yz,zx,yz,zx)- 
type orbital ordering. 

For the stabilization of this orbital ordering, the e g 
orbital degrees of freedom are indispensable. When the 
GdFe03-type distortion is absent, the hybridizations be- 
tween the neighboring t2 g and e g orbitals are prohib- 
ited by symmetry. As the GdFe03-type distortion in- 
creases, the symmetry restriction becomes to be relaxed, 
and the ti g -e g hybridizations are induced. Since the hy- 
bridizations between the neighboring e g and O 2p orbitals 
have a cr-bonding character and are large, the amplitudes 
rapidly increase by the distortion. Calculation of the 
spin-independent perturbational energy showed that the 
(yz,zx,yz,zx)-type orbital ordering is actually stabilized 
by the t2 g -e g hybridizations induced by the GdFe03-type 
distortion. Indeed, it was shown that a model without 
the e g orbitals cannot reproduce this orbital ordering. 

This orbital ordering was actuallyobserved in YTi03 
by several experiments (Refs. [H [H |U EJ IH 



indicating that this FM solution corresponds to nothing 
but the FM phase in YTi0 3 whose Ti-O-Ti bond angle is 
~ 140°. On the other hand, the AFM(A) phase has not 
been observed in the actual compounds so far. However, 
the nature of this FM phase is well characterized through 
the clarification of the AFM(A)-FM transition. Thus, we 
discuss the mechanism and nature of this AFM(A)-FM 
phase transition in the following. 

The AFM(A)-to-FM phase transition is identified as a 
transition where the sign of the spin-exchange constant 
along the c-axis changes from positive to negative while 
that in the afr-plane is always negative. The constant 
FM coupling in the a6-plane under the (yz,zx,yz,zx)- 
type orbital ordering can be easily understood. In the 
a6-plane, the neighboring orbitals are approximately or- 
thogonal to each other. Hence, the FM spin configura- 
tion is favored through the Hund's-rule coupling. How- 
ever, the emergence of the FM phase is not understood 
straightforwardly since the neighboring orbitals along the 
c-axis are not necessarily orthogonal. 

In fact, the hybridizations of the t?, g and e g orbitals 
again play an important role on the emergence of the 
FM phase. When the orbital state is strongly stabilized 
independently of the spin structure, the spin-exchange 
interaction is determined by the perturbational energy 
gain due to the electron transfers. In order to exam- 
ine the spin-exchange interaction along the c-axis, let us 
consider the perturbational processes with respect to the 
electron transfers between the 3d orbitals on sites 1 and 
3, which are neighboring along the c-axis. We fix the an- 
gles #afm(A) and #fm at 45°. Namely, we assume that 
an electron occupies -^(xy + yz), -^i~ x V + V z ) at sites 
1 and 3, respectively. 

In Fig. \7\ we display the main perturbational pro- 
cesses for both parallel and antiparallel spin configura- 
tions. When the spins are parallel, the transfer process 

2) is forbidden because of the Pauli principle. Therefore, 
the antiparallel spin configuration gains an energetical 
advantage of t\ju relative to the parallel spin configura- 
tion. On the other hand, in the transfer processes f ) and 

3) , the energies of the intermediate states are lowered 
by Hund's-rule coupling j. Therefore, the parallel spin 
configuration gains an advantage of, 



t 



1-2 



U'-J (U'-J+Aej 



+2 
l 3 



4i 



u' 2 + (u' + A e y 



(5) 



relative to the antiparallel configuration. Here, A e is 
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FIG. 8: Characteristic second-order perturbational energy 
gain due to the transfer processes along the c-axis for the 
antiparallel spin configuration t^/u, and that for the parallel 
spin configuration tfj/u' 2 + t$j/(u' + A Es ) 2 are plotted as 
functions of the Ti-O-Ti bond angle. Two energies are cross- 
ing at about 141°. With decreasing Ti-O-Ti bond angle, the 
sign of the spin-exchange constant along the c-axis changes 
from positive (AFM) to negative (FM) continuously at the 
phase boundary. (From Mochizuki and Imada, Ref. |34|'l 



the energy difference between the ti g and e g levels. 

Then, the spin configuration along the c-axis is de- 
termined by the competition between these two ener- 
gies, t\ju and t\j/u' 2 + t 2 j/{u' + A 6g ) 2 , which favor 
the FM and AFM spin configurations, respectively, and 
the spin-exchange constant is proportional to the differ- 
ence of these two energies. As the GdFe03-type distor- 
tion increases, the ti g -& g hybridization (£3) significantly 
increases while the hybridization between the neighbor- 
ing ^2g orbitals gradually decreases. Consequently, the 
energy t\j /u' 2 + t\j/(u' + A Cs ) 2 increases with increas- 
ing distortion while the energy t\ju decreases, leading 
to crossing of these two energies as shown in Fig |Sfa) . 
As a result, the spin-exchange along the c-axis changes 
from positive (AFM) to negative (FM) nearly continu- 
ously (see FigOJb)) We further note that since there ex- 
ists a very slight difference between the orbital states in 
the AFM(A) and FM phases, a tiny jump of the spin ex- 
changes should appear at the AFM(A)-FM phase bound- 
ary, resulting in the very weak first-order phase transi- 
tion. 

Now, let us summarize the above discussion. In 
the compounds with strong GdFeC-3-type distortion, the 
(yz,zx,yz,zx)-type orbital ordering is strongly stabilized 
irrespective of the magnetic structure due to both the d- 
type Jahn- Teller distortion and the t2 g -e g hybridizations 
induced by the GdFeOvtype distortion. Theoretically, 
if we assume this orbital ordering, the spin-exchange 
along the c-axis changes from AFM to FM with increas- 
ing GdFeOvtype distortion due to the increased t2 g -e g 
hybridizations, resulting in the AFM(A)-FM phase tran- 
sition. This indicates that for the emergence of the FM 
ordering in the perovskite titanates, the e g orbital de- 
grees of freedom play a decisive role even though the ti- 
tanates are the typical i2g-electron system. In addition, a 
nearly two-dimensional spin coupling with J c ~ is real- 
ized at the AFM(A)-FM transition point. Hence, in the 
FM phase, a two-dimensional anisotropy increases as the 
GdFe03-type distortion decreases. The decrease of Tq as 
the system approaches the AFM-FM phase boundary in 
the phase diagram is attributed to this two-dimensional 
anisotropy. Through these discussions, it is shown that 
the GdFe03-type distortion, which has been considered 
as a bandwidth-control mechanism so far, controls the 
orbital-spin structures in the perovskite titanates by in- 
troducing the t2 g -e g hybridizations. 

The importance of the e g orbitals can be easily un- 
derstood as follows. The order of the spin exchange due 



FIG. 9: Spin-wave spectrum for YTiOa measured by inelastic 
neutron-scattering experiment. The closed symbols refer to 
the magnetic excitations, while the open sym bols are identi- 
fied as phonons. (From Ulrich et aL, Ref. |4S|l 



to virtual transfer processes via the unoccupied e g or- 
bitals is t 2 /(u' + A Cg ) with t being a typical amplitude 
of the electron transfer between the Ti 3d orbitals. Here, 

the value of A e is estimated as ~ 1.7-2.0 eV from the 
1 1 1 1 1 1 1 I 

optical measurements k2l |43|, |4J, |45j and from the first- 
principles calculation [46(. On the other hand, that via 
the singly-occupied t2 g orbitals is t 2 /u = t 2 /(u' + 2j). 
Since A Cg is comparable to 2j ~ 1.3 eV, the e ff -orbital 
contribution t 2 / (v! + A es ) is not negligible, and actually 
plays a substantial role on the stability of the FM state 
in GdTi0 3 and YTi0 3 . 

Here, we note that the AFM ordering realized in the 
compounds with large R ions such as La, Ce, Pr, Nd and 
Sm, and that in the La-rich solid solutions Lai-^Y^TiOs 
is not A- type but G-type. In the above discussions, 
to focus on the FM compounds in the region of strong 
GdFe03-type distortion, the large d-type Jahn- Teller dis- 
tortion is assumed, which was observed in GdTiOs and 
YTi0 3 . Hence, the AFM(G) state in i?Ti0 3 (R =La, 
Pr, Nd and Sm) is beyond the scope of this framework, 
since in these compounds, the Jahn- Teller distortion is 
very small. 

The present calculations, which assume the strong d- 
type Jahn- Teller distortion show that the AFM (A) state 
is expected when the strong d-type Jahn- Teller distortion 
is realized in the region of moderate GdFe03-type dis- 
tortion. Here, "moderate" means "not strong enough to 
make the system FM" . However, this condition is, in fact, 
hard to be realized because the d-type Jahn- Teller distor- 
tion becomes weak rapidly with decreasing GdFe03-type 
distortion. Indeed, SmTi0 3 (ZTiOTi= 147°) does not 
show the d-type Jahn- Teller distortion while in GdTiO\3 
(ZTiOTi= 144°), a strong distortion is realized. There is 
no stoichiometric compound between SmTiOs (AFM(G), 
Ref. H3) and GdTi0 3 (FM) so that there is no possibil- 
ity of the AFM (A) ordering in the stoichiometric systems. 
However, this AFM(A) state may possibly be realized in 
the solid-solution systems like Lai-^Y^TiOs sandwiched 
by the AFM(G) and FM phases. The AFM(A) state may 
be found if we precisely control the magnitude of the 
GdFe03-type distortion by using the solid solutions or 
by applying pressures. However, it should be noted that 
even if the AFM(A) state is not realized in the actual 
systems, the validity of present discussion on the origin 
and nature of the FM state would never be affected. 

Finally, we mention a recent puzzling finding from 
the inelastic neutron-scattering experiment by Ulrich et 
al. (48J. They reported that the spin- wave spectrum of 
YTi03 is well described by the isotropic FM Heisenberg 
model (see Fig. |§J ■ On the other hand, if we take the 
(yz,zx,yz,zx)-type orbital ordering model, the FM cou- 
pling should have a two-dimensional anisotropy. This 
discrepancy has not been clearly resolved yet. To explain 
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the isotropic FM exchanges, Khaliullin and Okamoto 
recentl y p roposed a certain type of orbital ordering 
model [43, ■ This proposal is based on the analysis of 
a model with the triply degenerate ti g orbitals on the cu- 
bic lattice, and the effects of the actual lattice distortions 
(the GdFeC>3-type distortion and the Jahn- Teller distor- 
tion) are not taken into account. In fact, the predicted 
orbital structure contradicts the experimentally observed 
one. This contradiction may be attributed to the neglect 
of the lattice distortions, and may imply that it is neces- 
sary to take into account the lattice distortions when we 
consider the orbital-spin structures in the FM titanates 
as we have discussed in this section. This problem should 
be clarified by further studies in the future. Considering 
the fact that the energy scale of the phonon spectrum 
is comparable to that of the magnon spectrum, it may 
be necessary to consider the orbital excitations coupled 
with the lattice fluctuations. 



IV. AFM(G) ORDERING AND ORBITAL 
STRUCTURE IN LATI0 3 

LaTi03 with the smallest GdFe03-type distortion ex- 
hibits the AFM(G) ordering with T N -140-150 K. The 
magnetic ordered moment is 0.57 /is |12j . which is 
strongly reduced from spin- 1/2 moment of 1 /ib- Al- 
though the orbital-spin structure in LaTi03 has been 
studied intensively, even the origin of this AFM(G) or- 
dering had not been elucidated. However, a recent theo- 
retical study showed that the emergence of the AFM(G) 
ordering and puzzling experimental findings are con- 
sistently explained when we consider the crystal field 
from the La ions induced by the GdFe03-type distor- 
tion In this section, we would like to discuss the 
orbital-spin structure and several experimental results of 
LaTiOs on the basis of this theoretical proposal. 

Now, we first briefly summarize the history of stud- 
ies and controversies on the orbital-spin structure in 
LaTiOs. Previous diffraction studies showed that a Jahn- 
Teller type distortion of the Ti06 octahedron in LaTi03 
is undetectably small 0, E3- If this is true, the crys- 
tal field from the O ions surrounding the Ti 3+ ion has 
a cubic symmetry so that the ^-orbital degeneracy is 
expected to survive. 

Under this circumstance, the emergence of the 
AFM(G) ordering had been surprising since in the or- 
bitally degenerate systems, it is theoretically expected 
that a FM state with antiferro-orbital ordering is favored 
both by the electron transfers and by the Hund's-rule 

coupling mm m us m m . Indeed, a recent theo- 
retical study based on a Kugel-Khomskii model showed 
that in the cubic lattice with threefold degenerate ty. g 
orbitals, a FM state is almost always stable, and the 
AFM(G) state appears only in the unphysical regions of 
spin-exchange parameters [54j . 

On the basis of the Hartree-Fock analysis of the multi- 
band d-p model, Mizokawa and Fujimori proposed that 
the spin-orbit (LS) interaction lifts the ti g degeneracy, 



FIG. 10: Spin-wave dispersion in the (1,1,1) direction of the 
reciprocal space for LaTiOs measured by neutron-scattering 
experiment. The line is the fitting curve by using the Heisen- 
berg model with isotropic J of 15.5 meV on the cubic lattice. 
(From Keimer et al., Ref. [55f ) 



and the AFM(G) state in LaTiOs is accompanied by the 
LS ground state, in which two states -j^iv'z 1 + iz'x') | 

and ^(y' z ' ~ iz'x') J, with z'-axis and spins pointing in 
the (1,1,1) direction are alternating between the neigh- 
boring Ti sites 0, 0|. In the LS ground state, un- 
quenched orbital moment antiparallel to the spin mo- 
ment is necessarily induced. At first sight, this seemed to 
be consistent with the experimentally observed reduced 
magnetic moment in this compound. 

However, in contrast with this naive prediction, a re- 
cent neutron scattering study showed a spin-wave spec- 
trum well described by the isotropic spin-1/2 Heisenberg 
model (J ~15.5 meV) with a considerably small spin 
gap of —3 meV (see Fig. HU|) 55]. This indicates that 
the orbital moment is almost quenched, and the LS in- 
teraction is irrelevant since if the orbital moment exists, 
an induced Ising-type anisotropy in the spin sector has 
to generate a large spin gap. Actually, an exact diago- 
nalization study showed that the Kugel-Khomskii model 
with the LS interaction does not describe the observed 
spin- wave spectrum |56l . In addition, a model Hartree- 
Fock study done by Mochizuki j5]j also showed that the 
LS ground state cannot reproduce the AFM(G) ordering 
by examining a solution which was overlooked in the pre- 
vious Mizokawa and Fujimori's studies of Refs. [l4l il5| - 
In Ref. it was shown that even without the static 
Jahn- Teller distortion, a FM state out of which two states 
(yz + izx) 1 and xy f are alternating is stabilized both 
by the LS interaction and by the spin-orbital superex- 
change interaction. 

To explain the neutron scattering result, a possible 
orbital liquid state was proposed on the basis of small 
orbital-exchange interaction in the AFM(G) spin struc- 
ture HI, |H In these theories, the AFM(G) or- 
dering is assumed a priori. However, in the titanates, 
the spins and orbitals strongly couple each other, and 
both degrees of freedom cannot be determined indepen- 
dently. Therefore, the origin of the AFM(G) state in 
LaTiOs is to be clarified in a self-consistent manner. 
More importantly, a FM state with antiferro-orbital or- 
dering was theoretically expected to be more stable in 
this system. Thus, the assumption of the AFM(G) order 
on this basis is hard to be justified since relative sta- 
bility of the AFM(G) to FM state has never been ex- 
amined. In addition, it has been recently proven that 
the Kugel-Khomskii model with triply degenerate ti g or- 
bitals and zero Hund's-rule exchange on the undistorted 
cubic lattice, on which the orbital liquid theory is based 
does not exhibit the long-range magnetic ordering at any 
finite temperature |6lL 1(321 ] . This also questions a valid- 
ity of the orbital liquid theory since this theory assumes 
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FIG. 11: With the D^d distortion, the Ti06 octahedron is 
contracted along a trigonal direction. There exist two kinds 
of 0-0 bonds, shorter O-O bond and longer O-O bond as 
presented by solid lines and dashed lines, respectively. As a 
result of this distortion, the threefold degenerate t2 g levels 
split into a nondegenerate lower a\ g level and twofold degen- 
erate e g levels. With the D$d crystal field with [1,1,1] trigonal 
axis, the representation of the ai g orbital is -y=(xy + yz + zx). 



the AFM(G) ordering in the limit of weak Hund's-rule 
coupling. Indeed, a recent heat capacity measurement 
showed that the most of the low-temperature heat capac- 
ity arises from magnon contributions, which contradicts 
the prediction from the orbital liquid theory about the 
orbital contributions |r33| . 

A Raman scattering measurement performed for half 
metallic Lai_j,Ti03 with slight La-ion vacancies observed 
a Fano-type anomaly in the phonon spectrum [64|. The 
origin of this anomaly was identified as the coupling of 
phonons and orbital excitations in the orbital liquid the- 
ory 58, 59], and this observation was considered to sup- 
port this theory. However, another Raman scattering 
measurement performed for stoichiometric LaTiC>3 did 
not observe the anomaly |H3 • 

A possible trigonal (1)3,2) distortion of the TiOg oc- 
tahcdra was proposed to explain the emergence of the 
AFM(G) structure ^6]. With this distortion, the Ti0 6 
octahedron is contracted along the threefold direction, 
and the threefold degenerate ti g level splits into a non- 
degenerate lower a\ g level and twofold-degenerate higher 
e g level (see Fig. lllll . Occupations of the a\ g orbitals well 
explain the emergence of the AFM(G) ordering and the 
isotropic spin-wave spectrum. However, this distortion 
had not been observed clearly upon theoretical proposal. 

NMR measurements were also applied to this com- 
pound. Previously, the NMR spectrum was analyzed 
based on the orbital liquid picture [l^. However, it was 
recently claimed that there exists a discrepancy in this 
analysis. More specifically, the NMR spectra show the 
large quadrupole moment of a 3d electron in LaTiC>3, 
which is unfavorable for the orbital liquid model. In- 
stead, it was recently clarified that the spectrum is well 
described by an orbital ordering model proposed in the 
above Z^-distortion scenario instead of the orbital liquid 
model [67]]. Moreover, a recent resonant x-ray scattering 
result also indicates the orbital orderings in the series 
of the AFM(G) compounds of i?Ti0 3 (R =La, Pr, Nd 
and Sm), which have the same symmetry with the or- 
bital orderings observed in the FM compounds, GdTi03 
and YTi03 (68|. This also contradicts the orbital liquid 
picture. 

Several first-principles methods were also applied to 
LaTi03. However, LDA and LSDA failed to reproduce 
not only AFM(G) structure but also insulating behav- 
ior since the strong electron correlation is insufficiently 
treated HHHIzSo]. In addition, although the insulat- 
ing gap was reproduced, the AFM(G) structure was not 
obtained in the LDA+C/ calculation |70(. These indicate 



FIG. 12: (a) Displacements of the La cations around a TiOe 
octahedron (site 1) in the GdFe03 structure are presented by 
arrows. In LaO-plane 1, they shift in the negative direction 
along the b axis while in LaO-plane 2, they shift in the positive 
direction. The distances between the Ti ion and two La ions 
(gray circles) located in the ±(1, 1, Indirections are decreased. 
As a result, the La ions generate a crystal field similar to the 
Dz,d crystal field with [1,1,1] trigonal axis (see text). Note that 
rotation of the octahedron due to the GdFe03-type distortion 
is not presented, (b) Stacking of the TiOe octahedra as well 
as shifts of the La ions are presented. Note that Ti and O ions 
are not explicitly presented in this figure. (From Mochizuki 
and Imada, Ref. |51| 'l 



that in this compound, the effect of strong electron corre- 
lations is crucially important, which is of general interest 
of the condensed matter physics. 

In such ways as discussed above, the problem of the 
orbital-spin structure in LaTi03 had been a subject of 
hot debates, but a consistent theory had not emerged. In 
a number of previous studies, the degeneracy of the ti g 
orbitals had been assumed since the TiOg octahedra re- 
tain the cubic symmetry. However, this assumption is, in 
fact, a source of controversies. A recent theoretical study 
showed that the t2 g degeneracy is actually lifted by the 
crystal field from the displaced La ions in the GdFe03- 
type structure [5ll l52| . The orbital-spin structure and 
available experimental results can be well explained if we 
consider the effects of this La crystal field. 

Now, according to Refs. 0, |52] , let us discuss the ef- 
fects of the crystal field from the La ions, which are dis- 
placed by the GdFe03-type distortion in LaTi03. With 
the GdFeC>3-type distortion, the La ions shift mainly 
along the (l,l,0)-axis (namely, the &-axis), and slightly 
along the (1, — 1, 0)-axis (the a-axis) 0)Hd)lZ3l- There 
are two kinds of LaO-planes (plane 1 and plane 2) stack- 
ing alternatingly along the c-axis ((0, 0, l)-axis) as shown 
in Fig. ^] In plane 1, the La ions shift in the negative 
direction along the fr-axis while they shift in the positive 
direction in plane 2. Then, distances between Ti and 
La ions are no longer the same. For sites 1 and 2, the 
distances between the Ti ion and the La ions located in 
the ±(1, 1, Indirections decrease while those along the 
±(1, 1, — Indirections increase (see Fig. ^| (a)). On the 
other hand, the distances along ±(1, 1, — Indirections de- 
crease while those along ±(1, 1, l)-directions increase for 
sites 3 and 4. (In Fig. ^| (b), the stacking of each site 
is presented.) The changes of the other Ti-La distances 
are rather small. As a result, the crystal field from the 
La ions is distorted from a cubic symmetry. 

On the other hand, without any distortion of the TiOg 
octahedron, the crystal field from the ligand oxygens has 
a cubic symmetry. Under this circumstance, when we 
introduce the crystal field from the La ions, the threefold 
degeneracy of the cubic-i2g level splits. We expect that 
the La ions generate an attractive crystal-field potential 
and the Ti 3d orbitals directed along the shorter Ti-La 
bonds are lowered in energy because the nominal valence 
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FIG. 13: Orbital structure in the D^d crystal fields with 
[1,1,1], [1,1,1], [1,1,-1] and [1,1,-1] trigonal axes at sites 1, 2, 
3 and 4, respectively. The orbital wavefunctions at each site 

are ts^ xv + v z + zx )' 75 ( x v + v z + zx )> -^s( x v ~v z ~ zx ) 

and -j^(xy — yz — zx), respectively. Note that tiltings due 
to the GdFeC>3-type distortion are not presented explicitly in 
this figure, and the orbital basis xy, yz and zx are defined in 
terms of the x-, y- and z-axes attached to the tilted octahe- 
dron in reality. 

of the rare-earth ion is 3+. Therefore, we expect that 
the lowest orbitals are approximately ^( X V + U z + zx )i 
■^(xy + yz + zx), -^(xy-yz- zx) and ^= (xy - yz - zx) 
at sites 1,2,3 and 4, respectively. 

In order to confirm this, calculations based on the point 
charge model was performed by assuming the +3 valence 
on each La ion. In this calculation, the Coulomb inter- 
action between an electron on a Ti 3d orbital and sur- 
rounding La 3+ ions is given by using a dielectric constant 
£TiLa as, 

where Ri expresses the coordinate of the i-th La ion, 
which is deduced from the x-ray diffraction data Q , and 
Zr(=+3) is the nominal valence of the La 3+ ion. 

The calculation showed that due to the attractive 
Coulomb potential from the La 3+ ions, the ti g level splits 
into three isolated levels. The wavefunctions of the lowest 
orbitals at each site are specified by the linear combina- 
tions of the xy, yz and zx orbitals as follows, 

site 1; a\xy) + c\yz) + b\zx), 

site 2; a\xy) + b\yz) + c\zx), 

site 3; a\xy) — c\yz) — b\zx), 

site 4; a\xy) — b\y z) — c\zx) . (7) 

where a 2 +b 2 +c 2 =l. The coefficients take a=0.60, 6=0.39 
and c=0.69. Here, the xy, yz and zx orbitals are defined 
in terms of the x-, y-, and z-axes attached to each TiOg 
octahedra. The NMR spectrum proved to be well de- 
scribed by this orbital ordering model |73| . In addition, 
dynamical mean-field calculations combined with first 
principles local density approximation (LDA+DMFT) 
showed similar orbital ordering with nondegenerate or- 
bital ground state for LaTi0 3 [H Ell- 
in this orbital ordering, the orbital wavefunctions at 
sites 1 and 2 are actually similar to (xy + yz + zx) , 

and at sites 3 and 4 they are similar to A= (xy — yz — zx) . 
In fact, -^=(xy + yz + zx) is the lowest orbital for 
the D^d crystal field with [1,1,1] trigonal axis, while 
~^( X U ~ V z ~ zx) is the lowest orbital for that with 
[1,1,-1] trigonal axis. The orbital structure realized 
in the La crystal field is similar to that with the low- 
est orbitals in the D^d crystal fields with [1,1,1], [1,1,1], 



[1,1,-1] and [1,1,-1] trigonal axes at sites 1, 2, 3 and 4, 
respectively (see Fig. IT3)) . Therefore, this orbital order- 
ing is expected to induce the D%d distortions of the unit 
TiOg octahedra with these trigonal axes. In fact, the pos- 
sibility of the D$d distortions in LaTiC>3 was previously 
proposed in Ref. [66|. After this proposal, a number of 
diffraction studies were performed [12, [ZS Il2| > and such 
distortions were actually detected by an x-ray diffraction 
measurement |12|. 

One might argue that this TiC-6 distortion could occur 
by the usual Jahn- Teller mechanism. However, this is 
not true because the orbital degeneracy is already lifted 
by the La crystal field and the Jahn- Teller mechanism is 
not effective if the GdFe03-type distortion (tilting of the 
TiOe octahedra) is present. The observed TiC-6 distor- 
tion is simply induced by the orbital ordering due to the 
La crystal field. In this sense, the t^g degeneracy is lifted 
mainly by the La crystal field, and the Ti06 distortion 
occurs subsequently. 

Here, we note that the crystal field from the La ions 
has two origins. One is the Coulomb potential from the 
charged La ions as we have discussed above, and the other 
is the hybridizations between the Ti 3d orbitals and un- 
occupied orbitals on the La ions. The latter effect was 
also examined by using the second-order perturbational 
expansion in terms of the electron transfers between the 
Ti 3d and La 5d orbitals. The Slater-Koster parameters 
for the Ti 3d-La 5d hybridizations and the energy dif- 
ference between these two orbitals are deduced from the 
photoemission spectra and the LDA band calcula- 
tion |79| . respectively. It was shown that the hybridiza- 
tions between the Ti 3c? and La bd orbitals also stabilize 
the orbital ordering in Eq. l(7|l. and works cooperatively 
with the Coulomb attractive potential. 

We next discuss the stability of the magnetic state in 
the La crystal field (Hereafter, the Hamiltonian for the 
La crystal field is referred to as i?La)- An energy calcula- 
tion by using the effective spin-pseudospin Hamiltonian 
introduced in Sec. 3 was performed. Substituting the 
crystal field term H cry , with i?La, the energies for sev- 
eral magnetic structures were calculated by applying the 
mean-field approximation. This calculation showed that 
the AFM(G) ordering is the most stable solution. 

The spin-exchange constants were also calculated. 
Since the level splitting due to the La crystal field is 
much larger than fee^N, the orbital occupation is re- 
stricted to the lowest orbitals irrespective of the spin 
structure. Then, we can estimate the magnitude of 
the spin-exchange interaction in the subspace of singly- 
occupied lowest orbitals. The spin-exchange constant J 
for a Ti-Ti bond is represented as, J = (E^ — E^)/2S 2 
with and E^ being mean-field energy gains for the 
Ti-Ti bond of If- and ||-pairs, respectively. For LaTi03, 
the values of J along the x-, y- and z-axes take as 
,7^=18. 5 meV, 7^=18.5 meV and J z —19.7 meV, respec- 
tively, which well reproduce the spin- wave spectrum with 
isotropic exchange constant of ~15.5 meV. 

In addition, the ti g splitting is sufficiently larger than 
the coupling constant of the LS interaction in Ti 3+ 
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(Crf=0.018 eV) |s3 so that the La crystal field dominates 
over the LS interaction, resulting in the quenched orbital 
moment. This is also consistent with the observed small 
magnon gap in the spin-wave spectrum. 

As discussed above, by considering the i2 g -level split- 
ting due to the La crystal field, the emergence of the 
AFM(G) ordering and the spin- wave spectrum character- 
ized by (1) isotropic AFM coupling, (2) J ~ 15.5 meV, 
and (3)very small spin gap can be well explained. 

Here we mention that for more detailed description of 
the magnetic structure in LaTi03 , it is necessary to con- 
sider the Dzyaloshinskii-Moriya (DM) interaction. For 
example, the AFM(G) ordering in this compound is ac- 
companied by a weak FM moment due to spin cant- 
ing [81j, and this weak ferromagnetism is explained by 
the DM interaction [sif. In addition, the small magnon 
gap in the spin- wave spectrum is also attributed by this 
interaction |55|. However, the origin and overall features 
of the magnetic structure in LaTiC>3 are well explained 
even without this interaction. 

Now, let us discuss a problem of the reduced ordered 
moment in LaTi03. The ordered moment was measured 
as 0.57 hb in the neutron-diffraction study or as 0-45 
/ie hi the earlier powder neutron-diffraction study |8l| . 
These value had been considered to be much smaller than 
0.85 fiB of the spin-wave theory on the 3D Heisenberg 
model. This reduction was attributed to the antiparallel 
contribution of the orbital moment due to the LS inter- 
action. However, as discussed above, the recent neutron- 
scattering experiment showed that the orbital moment is 
fully quenched. Then, the origin of the moment reduc- 
tion has been recognized as a puzzle. However, this is not 
so puzzling actually when we consider the strong charge 
and orbital fluctuations in LaTiOs. 

In fact, LaTiC>3 has a small insulating gap of 0.2 eV, 
and this compound is located in the vicinity of the metal- 
insulator transition point 0, 01 ■ Thus, in LaTiC>3, a 
large amount of itinerant fluctuations of charges and or- 
bitals should remain. Although these fluctuations are 
completely neglected in the Heisenberg model in the insu- 
lating limit (in the limit of vanishing these fluctuations), 
we expect that in the actual compound, the fluctuations 
strongly reduce the magnetic moment. For example, a 
quantum Monte-Carlo simulation performed by White 
et al. showed that the charge fluctuations diminish the 
ordered moment of ~0.6 fiB for the 2D Heisenberg model 
to -0.35 ^ B for the 2D Hubbard model with U = At p|. 
In addition to the charge fluctuations, the imperfect nest- 
ing of the Fermi surface and the orbital fluctuations also 
work as reducing the magnetic moment [85l l86l l87l ] , and 
the observed reduced moment can be understood by con- 
sidering these effects as discussed in Ref. 52] in detail. 

We note that the insulating behavior and the magnetic 
properties such as magnetization and Tn for LaTiOs are 
highly sensitive to the slight off-stoichiometry induced 
by the excess oxygens or by the La-ion vacancies jg^, 
189 . l9(i| Indeed, the early measurement using a poorly 
stoichiometric sample observed the magnetic moment of 
0.45 hb HH, which is smaller than the recently reported 



value of 0.57 hb ^3- This difference may be attributed 
to the induced small amount of carriers, which cause the 
itinerant fluctuations. 

In this section, we have discussed controversies on the 
orbital-spin structure in LaTiOs. The La crystal field 
generated by the GdFeOvtype distortion lifts the Ti ti g - 
orbital degeneracy. This leads to the orbital ordering in 
the ground state instead of the orbital liquid, which pro- 
vides a consistent description for the emergence of the 
AFM(G) ordering and the puzzling experimental results. 
As we have discussed, the following experimental find- 
ings for LaTiOs are explained: (1) the spin- wave spec- 
trum with isotropic J of — 15.5 meV and small magnon 
gap measured by the neutron-scattering experiment, (2) 
the orbital structure detected by the NMR experiment, 
and (3) the trigonal distortions of the TiOg octahedra ob- 
served in the x-ray diffraction measurement. In addition, 
a quite recent spin-resolved photoelectron spectroscopic 
experiment with circular polarized light showed that in 
LaTiOs, the crystal field splitting of the t 2g subshell is 
about 0.12-0.30 eV and the orbital moment is strongly 
reduced js^. This is also consistent with the present 
theory. 

In the previous section, it has been shown that the 
GdFe03-type distortion has a universal relevance on the 
orbital-spin structures in the titanates through inducing 
the ti g -e g hybridizations. Now, we mention that the 
GdFe03-type distortion controls the orbital-spin struc- 
tures also through generating the R crystal field. It 
should be noted that this mechanism is not a usual Jahn- 
Teller one. In the Jahn- Teller systems, 3c? electrons on 
the transition-metal ions and the lattice degree of free- 
dom strongly couple each other, and the lattice spon- 
taneously distorts to lift the 3d degeneracy and lower 
the electronic energy. On the other hand, the GdFeOs- 
type distortion is caused by interactions between the R 
ions and O ions. More concretely, the octahedral tilt- 
ing is induced by energy gain in R-0 covalency bonding 
in the perovskite compounds with i?-site mismatching, 
and subsequent i?-site shifts occur by R-0 ionic interac- 
tions [9l|. This distortion exists irrespective of the 3d 
electronic state, and has nothing to do with any spon- 
taneous lift of the electronic degeneracy in contrast with 
the Jahn- Teller mechanism. The lifting of the degeneracy 
of the Ti 3d orbitals occurs just as a consequence of the 
crystal field from the La ions. In this sense, the mech- 
anism here proposed is different from the Jahn- Teller, 
which is inherent and universal with the GdFe03-type 
distortion. In fact, this mechanism also works in other 
AFM(G) titanates with R =Ce, Pr, Nd and Sm as we 
will discuss in the next section. 



V. ORBITAL-SPIN STRUCTURES IN i?TI0 3 
(R=LA, CE, PR, ND AND SM) 

In the previous section, we have discussed that the ef- 
fect of the La crystal field consistently explains the puz- 
zling experimental results as well as the emergence of 
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FIG. 14: Orbital wavefunctions at each Ti site are presented 
as functions of the Ti-O-Ti bond angle. The dotted line indi- 
cates 1 / VS. Open symbols at 140° represent the orbital wave- 
function of YTiOa obtained by polarized neutron-scattering 
experiments (a=V0A and c=VCU3, Ref. [UHi]). when the 
bond angle decreases, the orbitals continuously approach the 
character of YTi03. 



the AFM(G) ordering in LaTiOV In fact, it was shown 
that the crystal fields from the R ions also play impor- 
tant roles in i?TiC>3 with other R ions 0, • In this 
section, we discuss the orbital-spin structures in i?Ti03 
with R being La, Ce, Pr, Nd and Sm by considering the 
effects of the R crystal field. 

First, let us discuss the orbital state in the R crystal 
field of these compounds. The crystal field Hamiltonians 
based on the point charge model for these compounds 
were derived again by using the experimentally measured 
coordination parameters, and the wavefunctions of the 
lowest orbitals were calculated. The lowest orbitals for 
Hri of these compounds are expressed also by Eq. JTJl. 
In Fig. El we show the coefficients a, b and c as functions 
of the Ti-O-Ti bond angle. In LaTi03 with the smallest 
GdFe03-type distortion, the orbital wavefunctions are 
similar to those of the lowest orbitals in the D$d crys- 
tal fields, namely -^=(xy + yz + zx) or -^=(xy — yz — zx). 

The weights of the xy, yz and zx orbitals are nearly the 
same in LaTiC>3. 

As the GdFeOVtype distortion increases, the orbital 
wavefunctions become to differ from those of LaTiC>3 
since the locations of the R ions gradually deviate from 
the trigonal directions due to the increased tiltings. With 
increasing GdFe03-type distortion, the component of the 
zx orbital monotonically decreases at sites 1 and 3, while 
the yz-component decreases at sites 2 and 4. The re- 
ductions of these orbital occupations indicate that with 
increasing GdFe03-type distortion, the crystal field from 
the R cations tends to progressively stabilize the orbital 
structure in which sites 1, 2, 3 and 4 are approximately 
occupied by the orbitals as follows: 



site 1 

site 2 

site 3 

site 4 



-^( X V + zx ), 
— (xy-yz), 
-j=(xy- zx). 



In fact, this orbital structure is the same as the or- 
bital ordering which has been observed in YTi03 with 
large d-type Jahn- Teller distortion by several experi- 
ments Il8l E IH El, IH . More concretely, a and c 
seem to be extrapolated to values obtained by a polarized 
neutron study 0, E3 (a=V0A and c=V(h6). Further, 
a recent resonant x-ray scattering study showed that the 
orbital state continuously changes as R goes from La to 



FIG. 15: Spin-exchange constants along the x-, y- and z-axes 
are plotted as functions of the Ti-O-Ti bond angle. (From 
Mochizuki and Imada, Ref. 1521') 



Y, and both the orbital state in SmTi0 3 (AFM(G)) and 
that in GdTiOs (FM) have the same twofold symme- 
try, indicating that the symmetry of the orbital struc- 
ture does not change through the magnetic phase tran- 
sition |68j| . This experimental finding is consistent with 
the present theoretical result. 

Now, we mention a relation between the TiOg distor- 
tion and the R ions. In the previous section, we have 
mentioned that in LaTi03, the orbital ordering stabilized 
by the La crystal field favors the stacking of the trigonally 
distorted TiOg octahedra with [1,1,1] and [1,1,-1] axes 
alternatingly along the c-axis, which was actually de- 
tected experimentally. On the other hand, GdTiOs and 
YTi03 exhibit the d-type Jahn- Teller distortion in which 
the elongated axes of the octahedra are parallel in the c- 
plane as shown in Fig. EJ As a consequence of this distor- 
tion, the xy and yz orbitals are lowered in energy at sites 
1 and 3, while at sites 2 and 4, the xy and zx orbitals are 
lowered. In fact, the orbital wavefunctions in Eq. (JjJJ are 
represented by linear combinations of the twofold lower 
orbitals in the d-type Jahn- Teller distortion. This sug- 
gests that the d-type Jahn- Teller distortion in GdTiOs 
and YTi03 is favored by the orbital ordering due the 
R crystal field. On the other hand, a recent multiband 
d-p model study showed that with strong GdFe03-type 
distortion, the covalency between the O and R ions sta- 
bilizes the d-type Jahn- Teller distortion [j^- These two 
effects, namely, the orbital ordering due to the R crystal 
field and the R-0 covalency may cooperatively work to 
stabilize the d-type Jahn- Teller distortion. It is interest- 
ing to point out that in the titanates, the R ions in the 
GdFe03-type structure may control the distortion of the 
unit TiOg octahedron. 

We also mention that the orbital structure in SmTiOs 
may slightly deviate from the orbitals presented here be- 
cause of an additional effect of the O crystal field due to 
the relatively strong TiOg distortion. More concretely, 
we expect that in SmTiOs, the occupation of the xy 
orbital is somewhat increased at all Ti sites as will be 
discussed later. 

We next discuss the magnetic states in i?Ti03 with R 
being La, Ce, Pr, Nd and Sm. The AFM(G) spin struc- 
ture is also reproduced in these compounds. In Fig. 1151 
the calculated spin-exchange constants J x , J y and J z are 
plotted as functions of the Ti-O-Ti bond angle. This fig- 
ure shows that as the GdFe03-type distortion increases, 
the spin exchange gradually decreases, which is consis- 
tent with the gradual decrease of Tn as R goes from La 
to Sm. 

The Neel temperature Tn for SmTiOs is about 50 K. 
This is 40% smaller than Tn for LaTiOs. However, the 
calculated spin-exchange constant J x for SmTiOs is as 
large as 70% of that for LaTi03. This is puzzling when 
we consider the fact that Tn is proportional to J ^jj^ in 
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FIG. 16: Standard deviation from the mean for three inequiv- 
alent Ti-O bond lengths is plotted as a function of the Ti-O-Ti 
bond angle. The data of the Ti-O bond lengths are deduced 
from Ref. Q 



FIG. 17: (a) Energy-level structures of the crystal field Hamil- 
tonians Hl\ (upper panel) and Hl2 (lower panel) as func- 
tions of the Ti-O-Ti bond angle (see text), (b) TiO bonds 
in SmTi03 are elongated in the zy-plane, and are contracted 
along the z-axis. The data of the Ti-O bond lengths are taken 
from Ref. 0. (From Mochizuki and Imada, Ref. [5^|1 



the Heisenberg models. 

This discrepancy can be solved by considering the ad- 
ditional effect of the O crystal field. In the above anal- 
ysis, the R crystal field due to the GdFe03-type distor- 
tion is examined by assuming only the cubic O crystal 
field. However, the TiC>6 distortion in SmTi03 is rela- 
tively strong so that we expect a substantial contribution 
from the O crystal field to the orbital-spin structure. 

When the TiC>6 octahedron is strongly distorted from 
the cubic symmetry, the ligand oxygens are responsible 
for lifting of the orbital degeneracy through generating 
the crystal field. On the other hand, when the atomic 
TiOg distortion is small, it is expected that the R crystal 
field plays a substantial role on lifting of the t2 g and 
e g degeneracy, and dominantly determines the orbital- 
spin structures. Here, if the situation is inbetween, a 
competition between the R crystal field and the O crystal 
field occurs. 

The magnitude of the TiC>6 distortion is expressed by 
the standard deviation from the mean for the three in- 
equivalent Ti-O bond lengths (cr). In Fig. a obtained 
by the x-ray diffraction study |j| is plotted as a func- 
tion of the Ti-O-Ti bond angle. The value for PrTi0 3 
is deduced by the interpolation of LaTi03 and NdTi03. 
This figure shows that the values for i?=La, Pr and Nd 
are rather small, indicating relatively small Ti06 distor- 
tions. Thus, the R crystal field actually determines the 
orbital-spin structures in these compounds. On the other 
hand, the values for R=Gd and Y are rather large be- 
cause of the strong d-type Jahn -Teller distortions. As for 
SmTi03, the TiOg octahedra are moderately distorted, 
which leads to strong competition between the crystal 
field from the Sm ions and that from the O ions. In fact, 
the orbital-spin structure in SmTi03 has a two dimen- 
sional character due to this competition as discussed in 
the following. 

The t 2g level splittings due to the O crystal fields in 
i?Ti03 were estimated. The O crystal field Hamiltonians 
of the electrostatic origin (-Hli) were derived by consid- 
ering the Coulomb repulsive interaction between an elec- 
tron on the Ti ti g orbitals and the ligand 2 ~ ions within 
the point charge approximation. In addition, the O crys- 
tal field Hamiltonians for the hybridizations between the 
Ti t2 g and O 2p orbitals {H^) were also derived. Here, 
in order to study the effects of the TiO-bond length vari- 



ation, the O-Ti-0 bond angle is assumed to be 90°. In 
Fig-EK a )i the relative energies of each level for Hl\ (up- 
per panel) and those for Hl2 (lower panel) are plotted. 
Here, exio is a dielectric constant. In the compounds 
with small GdFe03-type distortion (i? =La, Pr and Nd), 
the differences of the energies for each level are rather 
small whereas in largely distorted GdTi03 and YTi03, 
the xy and yz orbitals are strongly lowered in energy at 
sites 1 and 3, and the xy and zx orbitals are lowered at 
sites 2 and 4 with large d-type Jahn- Teller distortion. 

In the case of moderately distorted SmTiOs, the TiO 
bonds are elongated along the x and y directions while 
along the z-direction, they are contracted in all TiOg 
octahedra. As a result, the O crystal field in SmTi03 
energetically lowers the xy orbitals, and consequently fa- 
vors the xy occupation at every Ti site, while both yz 
and zx orbitals are higher in energy. 

The increased a;y-occupation is expected to cause a 
two-dimensional anisotropy in the spin coupling in which 
the AFM coupling along the c-axis is rather weak. In- 
deed, the calculation of the spin-exchanges showed that 
when the strength of the O crystal field is comparable 
to, or is larger than the Sm crystal field, a nearly per- 
fect two-dimensional spin coupling is realized. The ob- 
served depression of Tn at SmTi0 3 may be attributed to 
this two-dimensional anisotropy. This two-dimensional 
anisotropy may be observed in SmTi03 . 

Now, let us summarize the present section. We have 
discussed that the R crystal field lifts the t2 g degeneracy 
also in i?Ti03 where R =La, Ce, Pr, Nd and Sm, and sta- 
bilizes the AFM(G) ordering as well as the orbital order- 
ing. The spin-exchange constant is gradually decreased 
with decreasing size of the R ion, which is consistent 
with the decrease of Tn in experiments. Further, with 
increasing GdFe03-type distortion, the orbital state in 
the R crystal field tends to continuously change towards 
the orbital state in the FM compounds such as GdTi03 
and YTi03 without changing its symmetry, which is in 
agreement with the resonant x-ray scattering results. In 
SmTiOs with moderate octahedral distortion, the Sm 
crystal field and the O crystal field compete strongly, 
which causes a two-dimensional anisotropy in the orbital- 
spin structure. The depression of Tn at R =Sm in the 
phase diagram is attributed to this anisotropy. With de- 
creasing size of the R ion, while the R crystal field dom- 
inantly determines the orbital-spin structure at R =La, 
Ce, Pr and Nd, the O crystal field becomes to domi- 
nate over the R crystal field at R =Gd and Y through 
the severe competition at R =Sm. As observed by sev- 
eral experiments, although the GdFe03-type distortion is 
rather large, the orbital structures in GdTi03 and YTi03 
with strong Jahn- Teller distortion are well characterized 
by the d-type Jahn- Teller distortion, which leads to a FM 
ground state. 

Here, we mention that for quantitative study on the 
orbital-spin structures in the perovskite titanates, it is 
necessary to treat the competition and cooperation of the 
crystal field from the R ions and that from the O ions in 
more precise way. However, because of the uncertainty 
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FIG. 18: Schematic magnetic phase diagram for i?Ti03 and 
solid-solution system Lai-^Y^TiOs in the plane of temper- 
ature and the GdFe03-type distortion. At R —La,, Ce, Pr 
and Nd with relatively small GdFe03-type distortions, the 
crystal field (CF) from the R ions dominantly determines the 
electronic structures, while in significantly distorted GdTiOs 
and YTiOs, the O crystal field dominates. In SmTi03, both 
crystal fields strongly compete. 



due to the dielectric constants exio and exm, it is difficult 
to estimate the contributions of these two crystal fields 
without ambiguity. This is left for the future studies. 

Through the above discussion, it has been clarified 
that in the perovskite-type transition-metal oxides, the 
R ions in the GdFe03-type structure generate the crystal 
field, which lifts the 3d orbital degeneracy. This mecha- 
nism of lifting the orbital degeneracy eventually plays an 
important role in controlling the orbital-spin structures 
in the perovskite compounds, competing with the usual 
Jahn- Teller mechanism. This mechanism has a general 
importance since the GdFe03-type distortion is a uni- 
versal phenomenon, which is seen in a large number of 
perovskite-type compounds. This mechanism may also 
play important roles on the electronic structures in other 
perovskite compounds. 



VI. PHASE DIAGRAM AND 
BANDWIDTH-CONTROL EXPERIMENTS FOR 
THE INSULATING TITANATES OTIO3 

We are now on the stage to discuss the whole nature of 
the exp erimentally obtained magnetic phase diagram 0, 
ITol IllL Il3| on the basis of the discussions in Sees. 2-4 
(see also Fig. H5)l. In this phase diagram, the coupling 
and interplay of the orbitals and spins are controlled by 
the GdFeC>3-type distortion, resulting in various orbital- 
spin phases and their phase transition as a function of 
this distortion. 

In LaTi03 with the smallest GdFeCvtype distortion, 
since the atomic distortion of the TiOg octahedron is 
small, the crystal field from the La ions dominantly lifts 
the tig degeneracy and determines the Ti 3d-orbital state. 
As a result, LaTiC-3 exhibits the AFM(G) ordering with 
the lowest orbital occupation in the nearly-trigonal La 
crystal field. The AFM(G) ordering due to the R crystal 
field is also realized in the compounds of Pr, Ce, Nd and 
Sm with relatively small octahedral distortion. In addi- 
tion, as R goes from La to Nd, Tn gradually decreases 
with decreasing spin-exchange constant. The decrease of 
the spin exchange is caused by the decrease of t 2g hy- 
bridizations via the O 2p states with increasing GdFeOv 
type distortion. 

In SmTiOs, the TiOg octahedra are moderately dis- 
torted so that the competition between the crystal field 
from the O ions and that from the Sm ions comes 
about. In SmTiOs, the O crystal field works as increas- 
ing the xy-orbital occupation in all Ti sites. As a re- 



FIG. 19: Experimentally obtained 2A act for _RTi03 is plotted 
as functions of U/W. Here, A act is the activation energy 
of the resistivity, and U /W is the strength of the electron 
correlation normalized to that of LaTiOs. The solid line is 
the result of least-square fitting with a linear function. (From 
Katsufuji, Taguchi and Tokura, Ref. |Tcj|') 

suit, the AFM(G) spin coupling becomes to have a two- 
dimensional anisotropy. The depression of Tn of SmTi03 
may be attributed to this anisotropy. 

In GdTiO-3 and YTi0 3 with strong GdFe0 3 -type dis- 
tortion, the electronic structure is dominantly deter- 
mined by the crystal field from the O ions with the largely 
distorted TiO-6 octahedra. These compounds show the 
FM ground state with an orbital ordering. This orbital 
ordering is stabilized both by the ci-type Jahn- Teller dis- 
tortion and by the ti g -& g hybridizations induced by the 
GdFe03-type distortion. The ti g -e g hybridizations also 
play an important role on the stability of the FM or- 
dering. We have shown that the FM coupling in the 
vicinity of the AFM-FM phase boundary has a strong 
two-dimensional character, resulting in the depressed To 

Near the AFM-FM phase boundary, both AFM and 
FM states have a two-dimensional anisotropy, which 
is expected to depress Tn and To The puzzling 
second-order like behavior of the AFM-FM transition 
can be understood from these depressions. In addi- 
tion, although there is no stoichiometric system between 
SmTi0 3 (AFM(G)) and GdTi0 3 (FM), the AFM(A) 
state might be observed in the solid-solution systems such 
as Lai-zYjTiOs and Gdi_. E Sm. E Ti03 sandwiched by the 
AFM(G) and FM phases. 

We have pointed out that the GdFe03-type distortion 
works as a control mechanism on the orbital-spin struc- 
tures in the perovskite titanates through introducing the 
t2 g -e g hybridizations or through generating the R crystal 
field. In addition, we have also pointed out that in the 
titanates, the t 2g degeneracy is inherently lifted by this R 
crystal field. This allows a single-band description of the 
low-energy properties of this system as a good starting 
point. 

Indeed, it has proved that the results of bandwidth- 
control experiments are well described by the theories 
based on the single-band Hubbard models. The band- 
width control can be achieved in the end compounds 
i?Ti03 or in the solid solutions like Lai-^Y^TiOs. Since 
the GdFe03-type distortion reduces the electron trans- 
fers between the neighboring Ti t 2g orbitals mediated by 
the O 2p states, we can control the Ti 3d bandwidth 
through this distortion by varying the R ions in i?Ti03, 
or in Lai-^Y^TiOs by substituting the La ions with 
Y ions. According to the tight-binding calculation, the 
bandwidth of the t 2g state decreases linearly with increas- 
ing magnitude of the GdFe03-type distortion, or with in- 
creasing Y concentration, and the bandwidth for YTiC\3 
is smaller by about 20 % relative to that for LaTiOs jg] • 

By controlling the bandwidth W, the evolution of 
the Mott-Hubbard gap as a function of the strength 
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FIG. 20: Calculated charge gap as functions of the on-site 
Coulomb repulsion U (a) for the single-band Hubbard model 
in the infinite dimension wi th se micircular density of states 
(From Georges et al., Ref. |l09|p . and (b) for that on the 
square lattice (From Watanabe and Imada, Ref. (§3)- ^ n 
Fig. (a), the solid line and the dotted line show the result 
obtained from exact diagonalization and that from iterated 
perturbation theory. In (b), finite-size charge gaps on L x L 
lattices are extrapolated to that in the thermodynamic limit 
given by the black square with a linear fit. The strength of the 
Coulomb repulsion U and the charge gap A are normalized by 
the bandwidth for (a), and by the nearest-neighbor transfer 
integral for (b). 

of electron correlation (U/W) has been investigated 0, 
U 0, |9^. The magnitudes of the gap were measured 
by extrapolating the optical conductivity spectra for 
Lai_2;Y x Ti03 0. The gaps for i?TiC>3 were also mea- 
sured from the optical conductivity and the temperature 
dependence of resistivity 0, El • It turned out that the 
gap A shows a linear increase as a function of U/W as 
shown in Fig. ^] The linear opening of the charge gap 
has been reproduced in several t heore tical studies on the 
single-band Hubbard models 0, Il09| . We show the cal- 
culated charge gap as a function of the electron interac- 
tion U for a single-band Hubbard model in the infinite 
dimension in Fig. 1201 (ah In addition, this linear behav- 
ior has been also reproduced for the single-band Hubbard 
model on the 2D square lattice as shown in Fig. [2U] (b) . 
These two results imply that the linear growth of the 
Mott gap is universal irrespective of the spatial dimen- 
sionality, if a nonzero transition point U — U c exists. The 
knowledge about the intrinsic lifting of the ti g degener- 
acy due to the GdFe03-type distortion in the titanates 
justifies the correspondence between the experiments and 
the theories based on the single-band Hubbard models. 

VII. FILLING-CONTROL METAL-INSULATOR 
TRANSITION 

As mentioned in Sec.l, the titanate system has been 
focused as a filling-control metal-insulator transition sys- 
tem. Hole doping to the Ti 3d bands is achieved in 
Ri- x A x TiOz+y by substituting the trivalent R ions with 
the divalent alkaline-earth ions (^4 =Ca, Sr, Ba), or 
by using the oxygen off-stoichiometry. By finely con- 
trolling the hole concentration in the systems such as 
Lai-^Sr^TiOa+y and Yi-^CazTiOa+j,, the nature of the 
transition has been intensively investigated, and the crit- 
ical behavior of the transition and the strongly filling- 
dependent properties near the transition point have been 
revealed. 

The nature of the filling-control metal-insulator tran- 
sition has been also theoretically studied. Most of these 
studies are based on the single-band Hubbard models al- 
though in reality, the titanate system with the Ti 3d or- 
bitals should be described by a multi-orbital Hubbard 
model. However, a lot of results obtained by the single- 



FIG. 21: Experimentally measured doping dependence of the 
specific-heat coefficient (7) and the magnetic susceptibility 
(x) in the metallic region of Lai-zSr^TiOa (lower panel), and 
the Wilson ratio xll normalized by 3/^|/7r 2 fc| (upper panel). 
(From Tokura et al., Ref. 3) 



FIG. 22: Calculated doping dependence of the specific-heat 
coefficient (7) obtained by the dynamical mean-field study 
on the infinite-dimensional Hubbard model. The experimen- 
tal result for Lai-^Sr^TiOs (Tokura et al., Ref. 3) is also 
presented fo r co mparison. (From Rozenberg, Kotliar and 
Zhang, Ref. |110| 'I Upper panel: Calculated Wilson ratio for 
the infinite-dimensional Hubbar d mo del as a function of dop- 
ing. (From Georges et al., Ref. \Wty ) 

band theories have shown good agreement with the ex- 
perimental results in spite of their extreme simplifica- 
tions. 

So far, it has been controversial whether the simple 
single-band description can be justified for this system, 
and a question why the experimental results are well ex- 
plained by the theories based on the single-band Hub- 
bard models has not been clearly answered. In the pre- 
vious section, we have mentioned that in the titanates, 
the R crystal field due to the GdFe03-type distortion 
eventually lifts the i2g-orbital degeneracy, which allows 
the single-band description of the ground-state and low- 
energy properties as a good starting point. Indeed, as 
shown in the previous section, the theories based on the 
single-band Hubbard models well reproduce the linear 
behavior of the charge gap growth as a function of U/W, 
which was experimentally observed both in the end com- 
pounds i?TiC>3 and in the solid solutions Lai-^Y^TiOs. 

In addition to these insulating systems, the inherent 
lifting of the t2 g degeneracy also allows the single- 
band description of the metallic hole-doped systems 
in the low-energy region. In this section, we review 
several filling-control experiments on the hole-doped 
systems Ri- x A x T\Oy, by comparing with the corre- 
sponding theoretical results. A measurement of the 
resistivity showed that the metal-insulator transition 
occurs at 1 ~ 0.05 for Lai-^Sr^TiOa 0, IH HI, at 
x ~ 0.2 for Nrb^TiOs [97|, and at a; ~ 0.4 for 
Y^^Ca^TiOs [IE Mill- The Hall coefficient R^ 1 
showed the existence of electron-type carriers (negative 
R^ 1 ), and estimated carrier density is proportional to 
filling n = 1 — x implying existence of a large Fermi 
surface 0, |9^|. Several measurements have shown the 
strongly filling dependent properties of the hole-doped 
metallic compounds near the metal-insulator phase 
boundary in (1) magnetic susceptibility (x) and specific 
heat (C = 7T), (2) resistivity (p), and (3) optical 
conductivity (<j(uj)). In the following, we discuss the 
measured results on the above quantities along with the 
theoretical results each by each. 

(1) Magnetic susceptibility (x) and specific heat (C = 
7 T) 
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FIG. 23: Temperature dependence of the resistivity (p) in 
metallic Lai-^Sr^TiOs. (From Tokura at al., Ref. Q) 

Tokura et al . measured the filli ng dependence of \ 
and 7 for La^^Sr^TiOa 0, llOOj . Both two quanti- 
ties are critically enhanced as the system approaches 
the metal-insulator boundary (x—0) from a metallic side 
with a nearly constant Wilson ratio Rw = x/l °f ^2 
against doping x (see Fig. I21[) . Similar e nhan cement 
was also observe d in Y^^Ca^TiOs |98l I9H llOOj and in 
Ndi-zCazTiOs |lf)lj |. These enhancements indicate a 
critical increase of the effective electron mass. Such a 
mass renormalization in Ri^ x A x TiO^, was also observed 
in several measu reme nts such as resistivity El , o ptical 
conductivity [98l Il0 2 | . Rama n scattering |9o. Il03l , pho- 
toemission spectra |104 Il05j . and NMR |l06j . In fact, 
the mass renormalization near the filling control tran- 
sition observed here is a generic feature of the filling- 
control metal-insulator transition, which was pre dicted 
theoretically by F urukawa and Imada (Ref. [107]) and 
Imada (Ref. [lOSf ). In addition, the experimentally ob- 
served enhancement of 7 is well reproduced by a dynami- 
cal me an-field st udy on the infinite-dimensional Hubbard 
model Il09t lllC| . In Fig. [221 the result of Rozenberg et 
al. (Ref. |ll0jp for 7 is depicted as a function of doping. 
In this calculation, to consider the Sr-doped La-titanatc, 
the value of U is chosen to be insulating but close to 
the metal-insulator phase boundary at zero doping. The 
experimental result of Tokura et al. for Lai-^Sr^TiOs 
(Ref. 4]) is also presented for comparison in the same 
figure. The figure shows a good agreement between them. 

As for the Wilson ratio, the calculated result shows 
roughly doping-independent behavior in the intermedi- 
ate doping region in agreement with the experimental 
observation. On the other hand, the calculated Wilson 
ratio vanishes at the transition point while experimen- 
tally it remains constant around 2. In addition, for 
a large doping region, the Wilson ratio becomes the 
nonintcracting value of ~1 in the calculation because the 
electron correlation is very weak in this region, which 
is also inconsistent with the experimentally observed 
constant behavior. These discrepancies may possibly 
be due to the limitation of the single-band description, 
and may be solved if we properly take into account the 
orbital degrees of freedom and the coupling between 
spins and orbitals. 

(2) Resistivity (p) 

Temperature dependence of the resistivity in 
Lai -^Sr^Ti Oa for various filling was mea- 
sured 0, IH 1 1 ll| . As shown in Fig. EH the resistivity 
is characterized by the relation p = po + AT 2 , indi- 
cating that the system is the Fermi liquid. This T 2 
behavior was r eproduced in the infinite-dimensional 
Hubbard model |l!2l Ill3| . Further, the coefficient A is 
enhanced as the system approaches the metal-insulator 
transition point [j. It was shown that a ratio A/7 2 
is approximately 1.0 x 10~ n fkm(molK/mJ) 2 , which 



FIG. 24: Experimentally measured optical conductivity spec- 
tra for Ei-xSr^TiCb and /Ji-^Ca^TiOs with R =La, Nd, Sm 
and Y. Each spectra is labeled by the nominal hole concentra- 
tion 8. Arrows show the equal-absorption points (see text). 
(From Katsufuji, Okimoto and Tokura, Ref. [j|) 



FIG. 25: The weight of the inner-gap spectra Nd as func- 
tions of the hole concentration S for .Ri-^Sr^TiOs and 
.Ri-zCazTiOs with R =La, Nd, Sm and Y. (From Katsufuji, 
Okimoto and Tokura, Ref. |Sj) 

is nearly the same as the universal value for the heavy 
fermion systems. This value is in fair agreement with the 
theoretically pre dicte d value in the infinite-dimensional 
Hubbard model |114| . 

(3) Optical conductivity (<r(u))) 

The evolution of the Drude part of the optical conductiv- 
ity as a function of doping also shows agreement between 
theories and experiments. The Drude weight is defined 
as the singularity of the conductivity at u> —0 and T =0 
in the strict sense. However, in reality, the experimental 
results always suffer from the finite-temperature broad- 
ening as well as life-time effects from the impurity scat- 
tering. Thus, experimentally the Drude weight is often 
discussed based on the inner-gap spectral weight of the 
conductivity. 

The doping-induced changes in the optical conductiv- 
ity for i?i_ iC Sr ie Ti03+y and Ri^ x Ga. x Ti0^ y were inves- 
tigated H, Hi The optical conductivity for R =La, 
Nd, Sm and Y are displayed in Fig. |^] The spectra are 
labeled with the nominal hole concentration 8 = x + 2y 
where y is the oxygen off-stoichiometry. These fig- 
ures show that the spectral weights transfer from the 
Mott-Hubbard-gap region to the inner-gap region with 
hole doping across an equal-absorption (isosbetic) point 
around 1.2 eV, leading to the growth of the quasi-Drudc 
weight. The effective number of electrons defined as 

N eS (u) = ^- f ' o{J)dJ (9) 

where mo is the free electron mass and N is the number 
of Ti atoms per unit volume, is nearly independent of 6 
with u> =2.5 eV, which includes the spectral weight of 
both the Mott-Hubbard gap excitations and the inner- 
gap excitations. This implies that the sum of these two 
parts is nearly conserved irrespective of 5. Therefore, the 
quantity N D = N e g(uj c )s ~ Acff.(w c )<5=o with u c ~1.1 eV 
being the isosbetic point denotes the lower-energy spec- 
tral weight transferred from the higher-energy region, or 
in other words, quasi-Drude weight. In Fig-ESl the values 
of No for R =La, Nd, Sm and Y are plotted as functions 
of hole concentration 6, which show the linear increase. 

On the other hand, the optical response was theoreti- 
cally calculated based on the single-band Hubbard model 
in the infinite dimension by using the quantum Montc- 
Carlo method combined with the dynamical mean-field 
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FIG. 26: Doping dependence of the optical conductivity in 
the infinite-dimensional Hubbard model with U — At on the 
hypercubic latti ce w ith t/(2y/d) from Jarrell, Freericks and 
Pruschke (Ref. | 115j| 1. The curves correspond to S =0.068, 
0.0928, 0.1358, 0.1878, 0.2455. 0.3, 0.4 and 0.5 from top to 
bottom curve at high u). The inset shows the evolution of the 
Drude weight as a function of doping, which shows a linear 
increase with doping. 



FIG. 27: Total spectral weight of the optical conductivity as a 
function of the electron concentration n = 1 — 5 in the infinite- 
dimensional Hubbard model for the same parameters as in 
Fig. 1261 The weight is decomposed into a Drude weight part 
(open triangles), a mid-infrared part (open squares) and a 
charge-transfer part (sol id triangles). (From Jarrell, Freericks 
and Pruschke, Ref. (lT^.I 

theory |115|| . We show in Fig. [23 the optical conduc- 
tivity sp ectra for small dopings obtained by Jarrell et 
al. (Ref. |ll5| ). They decomposed the spectral weights 
into the following three different contributions; (1) Drude 
part from the transitions within the quasiparticle reso- 
nance, (2) mid-infrared part from the transitions from 
the lower Hubbard band to the unoccupied part of the 
quasiparticle peak, and (3) charge transfer part from the 
excitations between the lower and upper Hubbard bands. 
The weights are presented in Fig. |57| where the Drude 
part linearly increases with increasing hole number 8 
(decreasing electron concentration n = 1 — 5), and the 
weights associated with the incoherent contributions are 
concomitantly decreased. The calculation shows a rea- 
sonable agreement with the above experimental results 
of Ref. @. 

As shown in this section, the single-band theories have 
succeeded in describing a number of results of the filling- 
control experiments. The intrinsic splitting of the ti g 
levels due to the R crystal field justifies the single-band 
description and provides a firm ground for the correspon- 
dence between the experiments and the single-band theo- 
ries. We expect that Lai-xSr^TiOa+y as well as LaTiC>3 
are especially well described by the Hubbard model on 
the cubic lattice where amplitudes of the electron trans- 
fers along the x-, y-, and z-directions are nearly the 
same since in these compounds, the GdFeC>3-type dis- 
tortion is relatively weak, and hence the lowest orbitals 
are directed nearly along the trigonal directions, which 
generates nearly isotropic electron transfers between the 
neighboring lowest orbitals. 

However, in YTiC>3 and GdTiC>3, the lifting of the or- 
bital degeneracy is strongly coupled to the Jahn- Teller 
distortion in contrast to the case of LaTi03 where the 
Jahn- Teller coupling is weak, and is dominated by the 
R crystal field. In this case, the orbital order may 
be strongly modified or even lost when the carriers are 
doped. Such complexity with contribution from the or- 
bital degrees of freedom combined with dynamical Jahn- 
Teller fluctuations has to be treated separately beyond 
the single-band scheme. In fact, the insulating behavior 



is extended to a heavily doped region in YTiC>3, where 
the orbital order and Jahn- Teller distortion play crucial 
roles. 

Finally, even for LaTi03, it may be necessary to take 
into account the orbital degrees of freedom for higher- 
energy excitations. For instance, the multi-orbital ef- 
fect may appear in the higher-energy region of the op- 
tical conductivity and the photoemission spectra. In- 
deed, it was recently claimed that the coherent and in- 
coherent features of the observed photoemission spec- 
tra for the doped titanates show disagreement with the 
results obtained from the single-ba nd Hubbard mod- 
els [H EH1 [To! EH EH Ell EH • To solve the dis- 
crepancies, these experimental data should be discussed 
in the light o f theories based on the multi-orbital Hub- 
bard models |l2Ct Il2 1| . The multi-orbital effects in the 
higher-energy excitations are issues of interest, and fur- 
ther studies are needed. 



VIII. CONCLUDING REMARKS 

The perovskite Ti oxides have been recognized as a 
key class of materials for understanding the physics of 
strongly correlated electrons because of the 3d 1 configu- 
ration of the end compounds. The coupling of magnetism 
and orbitals is one of the most interesting subjects to be 
clarified in this field, and the titanates have provided an 
important example for this subject. A rich magnetic and 
orbital phase diagram has attracted great interest. 

In this article, we have reviewed both theoretical and 
experimental studies and current understanding on the 
orbital-spin structures in the perovskite titanates as well 
as their phase transitions. Through these discussions, we 
have shown that there are several control mechanisms, 
which determines the orbital-spin structures in the ti- 
tanates. The obtained knowledges about the orbital-spin 
couplings and the control mechanisms are expected to be 
helpful for clarification of electronic structures in other 
perovskite compounds. For instance, we have discussed 
that the GdFe03-type distortion has a universal rele- 
vance on the orbital-spin structures through introducing 
the ti g -e g hybridization (cf. Sec. 3) or through gen- 
erating the R crystal field (cf. Sees. 4 and 5). The 
mechanisms here proposed may work in other perovskite 
compounds since the GdFe03-type distortion is a phe- 
nomenon which is seen in a large number of perovskite 
compounds. Thus, we expect that several puzzling prop- 
erties and phenomena of other perovskite compounds can 
also be clarified in the future research if we consider these 
mechanisms. The knowledges may be important also 
from the aspect of material design. Manipulation of the 
coupling of spin, orbital and lattice degrees of freedom 
through these control mechanisms may open a possibil- 
ity of controlling the dimensionality, transport properties 
and cooperative response to the external fields of these 
materials. 

In addition to the coupling of magnetism and orbitals, 
the metal-insulator transition (Mott transition) is an- 
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other important subject of the physics of strongly corre- 
lated electrons. There are two different routes to achieve 
the Mott transition. One is the bandwidth control (con- 
trol of the electron correlation strength U/t) and the 
other is the filling control. As we have discussed in this 
article, the titanate system provides typical examples of 
these issues as well. The filling-control Mott transition 
in the hole-doped titanates has been intensively stud- 
ied. On the other hand, the bandwidth control can also 
be achieved in the insulating titanates, and the charge 
gap evolution as a function of bandwidth has been in- 
vestigated although the bandwidth cannot be sufficiently 
large to make the system metallic. 

We have discussed the results of bandwidth-control ex- 
periments and those of filling-control experiments in the 
light of the theoretical results based on the single-band 
Hubbard models. The single-band description of this 
system as a starting point is justified as far as for the 
ground-state and the low-energy properties since the Ti 
i2 S -orbital degeneracy is inherently lifted by the R crys- 
tal field in the GdFe0 3 -type structure. It has turned out 
that the titanate system embodies the essential physics 
of the single-band Hubbard model which is one of the 
simplest and the most important model capturing the 
competition of the electron interactions and the kinetic 



energies in the strongly correlated electron systems. 

As discussed in this article, the intensive studies have 
elucidated the ground-state electronic structures of the 
perovskite titanates. On the other hand, clarification of 
the properties of excited states is one of the important 
subject of further works. For instance, the photoemission 
spectroscopy and the optical conductivity measurement 
are expected to detect the higher-energy excitations. To 
analyze the spectra obtained in the experiments, it may 
be necessary to take into account the multi-orbital ef- 
fects on the excitations. Indeed, the filling dependence 
of photoemission spectra shows a discrepancy between 
the experiments and the single-band theories as we have 
mentioned in Sec. 7. It is interesting to study the spec- 
tra based on the multi-orbital Hubbard model. In ad- 
dition, the possibility of detecting the collective orbital 
excitation (so-called "orbital wave" ) in th e titanates has 
recently been proposed theoretically 122]. Prediction of 
the nature of orbital-wave spectrum or its experimental 
detection are also interesting. For the future research on 
such properties and phenomena associated with the exci- 
tations, the understanding of the ground-state electronic 
structures discussed in this article will be essentially im- 
portant. 
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